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Incompressible polaritons in a flat band 
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We study the interplay of geometric frustration and interactions in a non-equilibrium photonic lattice system 
exhibiting a polariton flat band as described by a variant of the Jaynes-Cummings-Hubbard model. We show 
how to engineer strong photonic correlations in such a driven, dissipative system by quenching the kinetic energy 
through frustration. This produces an incompressible state of photons characterized by short-ranged crystalline 
order with period doubling. The latter manifests itself in strong spatial correlations, i.e., on-site and nearest- 
neighbor anti-bunching combined with extended density-wave oscillations at larger distances. We propose a 
state-of-the-art circuit QED realization of our system, which is tunable in-situ. 


Over the last decade there has been a surge of interest in 
realizing strongly correlated states of light in interacting pho¬ 
tonic lattices for quantum simulations and the study of non¬ 
equilibrium many-body physics QHEl (for two recent re¬ 
views, see [13! [f4|). Effective photon-photon interactions 
can be engineered in these systems utilizing strong light- 
matter couplings in various cavity/circuit QED platforms, e.g., 
with atoms lfl5l . excitons 1161 [T7j | or superconducting qubits 
mm. Arranging cavities and atoms/qubits on a lattice of¬ 
fers the opportunity to engineer strongly correlated states of 
photons in various geometries with local control over coher¬ 
ent as well as dissipative dynamics. The driven dissipative 
nature of photonic systems then allows for direct experimen¬ 
tal, non-invasive access to the complete density matrix, e.g., 
temporal and spatial correlation functions l20l . 

A particularly challenging and interesting problem of 
many-body physics concerns the study of frustrated lattices. 
Frustration refers to the impossibility of satisfying simulta¬ 
neously all constraints implied by a Hamiltonian, which are 
imposed, e.g., by geometry, disorder or interactions. This 
typically leads to macroscopically degenerate ground-states, 
which are sensitive to small perturbations and thus define a 
challenging minimization problem. Conversely, frustration 
often gives rise to interesting strongly correlated phenomena 
and the emergence of fascinating non-trivial structures, e.g., 
in quantum magnetism BTl - BSl . quantum hall systems 126! 
|28l , Josephson junctions 129 301 or ultra-cold atoms I13lti 35lh 
In this work, we make use of geometric frustration to boost 
interactions and show that photons pumped into the flat band 
of a photonic lattice form an incompressible state of light with 
non-trivial spatial correlations at the onset of crystallization. 
This steady state cannot follow from energy minimization, but 
originates under non-equilibrium conditions with balanced 
drive and dissipation. 

First realizations of interacting photonic lattices have re¬ 
cently been engineered based on superconductor as well as 
semiconductor technologies Ii36fl40l . Motivated by these 
achievements, we study a ID qubit-cavity chain, where qubits 
couple to photons in every other cavity. Such a Jaynes- 
Cummings-Hubbard (JCH) system |2] [3] [4lfl50l can be 
viewed as a quasi-ID cut through a 2D Lieb lattice GU, 
where the qubit represents one of the sites in the unit cell 
and simultaneously generates the frustration leading to the flat 
band as well as the photon-photon interaction. This setup is 
readily realizable with state-of-the-art circuit QED technol¬ 


ogy, where the lattice dispersion as well as the strength of the 
effective photon-photon interactions can be tuned in-situ by 
simply changing the qubit-resonator detuning DU. We show 
that one of the Bloch bands of the array can be tuned from 
dispersive to completely dispersionless, i.e., flat in the entire 
Brillouin zone. This flat band arises due to destructive quan¬ 
tum interference and generates a macroscopic set of degener¬ 
ate and localized plaquette states Il22ll23i (further information 
on creating flat bands can be found in the supplemental mate¬ 
rial (SM), which includes Refs. EH1). Similar flat bands 
were recently observed in non-interacting 2D laser and micro¬ 
pillar arrays EBlISTl. 

Here, we investigate the effects of strong photonic interac¬ 
tions in such a non-equilibrium geometrically-frustrated sys¬ 
tem by using projective methods as well as a time-evolving 
block decimation (iTEBD) algorithm l58l [59 1. We find that 
geometric frustration strongly enhances photon repulsion on 
the lattice and pushes the system towards an incompressible 
state characterized by short-ranged crystalline order with pe¬ 
riod doubling. Incompressibility is signaled by the appear¬ 
ance of an extended plateau in the average polariton excitation 
number as a function of drive strength, whose height is deter¬ 
mined solely by the geometry of the lattice. Crystallization 
manifests itself in strong spatial photonic correlations, i.e., 
on-site and nearest-neighbor anti-bunching combined with ex¬ 
tended density-wave oscillations at larger distances. Interest¬ 
ingly, we find that the correlation length of these oscillations 
can be increased when decreasing the light-matter coupling 
strength g with respect to the photon hopping rate J. 



FIG. 1. (color online). Top: Sketch of a transmission line resonator 
chain including superconducting qubits in every other resonator. Bot¬ 
tom: Simplified lattice representation. Two resonators of type A and 
B (rectangles) are coupled by the photon hopping rate J. Qubits (cir¬ 
cles) are only coupled to the A cavities with strength g. The dashed 
lines show one unit cell of the array and the blue symbols mark a 
localized plaquette state as discussed in the text (± denote the corre¬ 
sponding phases in the wavefunction), see Eq. 0 - 
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We study a variant of the driven dissipative JCH model, i.e.. 
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where hj a denote the on-site Hamiltonians for resonators of 
type A and B with qubits at site Q coupling only to the A 
sites, i.e., h jA = A A o]aj + A q <t+ a~ + (ga^aj+faj +H.c.) 
and hj B = A B bjbj + f(bj + H.c.). The bosonic operators 
cij ( bj ) annihilate a cavity photon at site A (B) in unit cell 
j = 1 ,,N. The second term in (0 describes photon hop¬ 
ping between nearest neighbor resonators at a rate J. The 
qubits are represented by Pauli operators crj and couple to 
the A photons with strength g. All cavities are subject to a co¬ 
herent drive of strength / described by the last terms in hj a . 
In a frame rotating with the drive frequency w D the bare cavity 
and qubit frequencies w x are renormalized to A x = — w D , 

with X = A, B, Q. Cavity dissipation is taken into account 
using a Lindblad master equation for the density matrix, i.e., 
p = —i[H, p] + (k/2) J2j(P[ a j]P + V[bj]p), with the Lind¬ 
blad operator V[a\p = 2 apa) — a^ap — pa 11 a and the photon 
decay rate n. Here we neglect spontaneous emission and de¬ 
phasing of the qubits, which can be substantially suppressed 
with respect to cavity decay lf60l . Fig. [Tjshows an implemen¬ 
tation of our model based on state-of-the-art circuit QED tech¬ 
nology m. A similar geometry can be realized using semi¬ 
conductor micro-pillar arrays, where the qubit site is replaced 
with a nonlinear cavity (see the SM fur further information) 

mm. 

We start with the discussion of the single-particle spectrum 
of {1} in the absence of drive (/ = 0) and dissipation (k = 0). 
For that purpose, we write a common Fourier transform 1 /> ■ = 
(1/VN) Ylk with i/jj = [dj, bj, a~ ] T and impose 

periodic boundary conditions to obtain the fc-space represen¬ 
tation of the lattice Hamiltonian, i.e., H = Y2k 
with k = 2nir/N, n = —N/2,...,N/2 — 1. The eigen¬ 
value equation for hk yields three bands, which are plotted in 
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FIG. 2. (color online), (a) Single particle dispersion of the lat¬ 
tice in Fig. [T] For <5q B = J (dashed) all three bands are disper¬ 
sive. For <5q B = 0 (solid) the middle band (blue) becomes flat 
corresponding to a set of degenerate single-particle plaquette states, 
(b) Many-body eigenstates associated with the flat band are con¬ 
structed from products of plaquettes (see main text). They form an 
equally spaced multi-level system indexed by the particle number 
n — 0,.., n m a X with degeneracies d n = ( 2rlnl " _71 ) anc j ingest filling 
n max = (N + l)/2. Here we show the case for N = 13 unit cells. 


Fig. |2ja). For the general case with 5 QB = w Q — w B ^ 0 
all bands are dispersive (dashed lines). However, if 5 QB = 0 
the middle band turns flat with energy e PB = w B while the 
other two remain dispersive with energies e± = w B +5 AS /2± 
y / 2J 2 ( 1 + cosfc) + g 2 + 6 AB /4, where S AB = w A — w B . The 
flat band eigenstates can be written as 

I A j> = y g2 ^ 2J2 [9 b ] - J K + + <+i)] l vac >» ( 2 ) 

which describes a localized plaquette state defined by one B 
and two neighboring Q sites (see Fig. [T]t. The flat band arises 
due to the destructive interference between a photon hopping 
process from resonator B to A (~ J) and the excitation trans¬ 
fer due to the coupling of qubit Q to the resonator A (~ g). 
As a consequence the A cavities remain completely dark, such 
that an excitation originally localized at one end of the chain 
does not disperse and/or propagate to the other end. 

In the following, we are interested in the interplay of frus¬ 
tration and interactions in the non-equilibrium steady state 
(NESS) of the system. In equilibrium, thermalization would 
lead to a zero temperature ground state not involving flat band 
states as these do not reside at the lowest energy. In a non¬ 
equilibrium setup however, the coherent drive can excite the 
flat band by keeping the drive frequency resonant with the flat 
band energy, i.e., oj d = e FB . In order to take into account 
states resonant with the drive, we construct from Q the eigen¬ 
states of the Hamiltonian {!]) (for vanishing drive amplitude, 
i.e., / = 0) with energies that are integer multiples of e FB , 
and project the Lindblad equation on this subspace. Apart 
from the single-particle states in |2|, these states are prod¬ 
ucts of non-overlapping plaquettes, e.g., the two particle states 
1 ^ 2 ) ~ |Ai) |A 3 ), |Ai) |A 4 )... with energy 2e FB , the three 
particle states | ip 3 ) ~ |Ai) |A 3 ) |A 5 ),... with energy 3e FB 
etc. (for further details see the SM). The energetically high¬ 
est lying state is the density-wave |d/dw) = n”:T 1 a 2j — 1 ) 
with energy £d w = n ma x£F B an d particle number n max = 
(N + l)/2, i.e., filling per lattice site v j w = n max /(3N) = 
1/6 + 0(1/N). This special ladder of flat band states with 
degeneracies d n = ( 2 ” m “ -n ), where n is the particle number 
of each state, is shown in Fig. [2jb) for N = 13 unit cells. 
All eigenstates with n > n TTmx belong to dispersive bands and 
are gapped from the flat-band ladder due to the nonlinearity 
induced by the light-matter coupling g. 

Due to the coherent drive with w D = e FB we expect states 
belonging to the flat band ladder to mostly contribute to the 
NESS at small and intermediate drive strength. In Fig. [3] we 
show the average excitation number per lattice site (v) (for a 
formal definition see caption of Fig. [3) as a function of pump 
strength f /n. At weak pump / <C k, the results of the pro¬ 
jected model (solid blue line) agree with the analytical expres¬ 
sion (v) « (4/ 2 /k 2 )[1 + (4J 2 + K 2 /A)/g 2 ] (straight solid 
line), which is obtained from a perturbative calculation of the 
steady state to leading order in // n. At stronger pump, how¬ 
ever, the system saturates at a filling (y) « r , d w /2 ~ 1/12 
resulting in an extended plateau centered around / ~ k. This 
plateau can be interpreted as an incompressible state of pho¬ 
tons with d(y)/df ~ 0, as we now explain in more detail. 
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FIG. 3. (color online), (a) Excitation number (v) = Ex( z/x )’ with 
vx = n x / (37V), n x = artjcj- (X = A, B. Q and Xj = a,j,bj, a~) 
in the steady state as a function of pump strength //ft. Shown are 
results obtained from projection of the density matrix on the flat band 
eigenspace for a system with IV = 13 unit cells and open boundary 
conditions (solid line) and from iTEBD simulations of the infinite 
system at zero detuning 5 q B = 0 (blue symbols) and finite detuning 
(5 qb = J (circles). The plateau at (i/) « v ~ 1/12 is associated 
with a suppression of number fluctuations K = [((Ex^) 2 ) — 
(Ex( nx )) 2 ]/Ex( nx ) as shown in (b). The plateau is extended 
for the difference {v) — (i / A } (asterisks), but almost vanishes in the 
dispersive case (circles), (c) Probability p„ of finding n excitations 
in the lattice as calculated within the flat band model for the pump 
strengths marked with arrows in (a). Other parameters: g/J = 1, 
5ab/J = 0.5, k/J = 0.05. 


The height of the plateau is largely independent of g and J 
and determined mostly by the geometry of the lattice. This 
can be understood by looking at the excitation number dis¬ 
tribution p n of finding n excitations in the lattice, shown in 
Fig. [3jc). At weak pumping the distribution is peaked at 
low excitation numbers and shifts to larger n for increasing 
pump strength. At strong pumping it saturates and resembles 
approximately the degeneracies d„ shown in Fig. [Tj'b), i.e., 
all states are almost equally occupied similar to a two-level 
system saturating half way between ground and excited state 
EEL The saturated average excitation number is thus calcu¬ 
lated as n « (E^o«d„)/(E^ 0 d„) = (1 - l/V5)(N/2) 
corresponding to roughly half the density-wave filling v = 
n/(3N) « r'<i w /2 (horizontal dashed line in Fig. |3ja)). The 
incompressible state thus originates from an unconventional 
photon blockade on a frustrated lattice arising from a satura¬ 
tion of the flat band ladder shown in Fig.[2jb). 

We confirm this picture by numerical simulations employ¬ 
ing an open system version of the iTEBD algorithm ll58l l59l 
(for technical details see SM). In Fig. |3|a) the projected 
model agrees with the exact numerics (diamonds) well into 
the plateau, thus verifying the incompressible state of pho¬ 
tons, where fluctuations of the excitation number are reduced 
(see Fig. |3jb)). For an even stronger pump (/ ft), the dis¬ 
persive bands start to contribute to the NESS leading to a de- 



FIG. 4. (color online). Correlation function of photons emitted by the 
B sites = (bQbjbobj) /(b^bo) (bjbj) for different drive strength’s 
//k at fixed detuning <5q B = 0 (upper panel) and different detunings 
<5 qb/ J at fixed drive strength //k = 0.05 (lower panel) as calculated 
with iTEBD. The density-wave oscillations correspond to a period 
doubling with respect to the unit cell of the underlying lattice. In (b) 
the drive stays resonant with the top of the middle band. The inset 
shows the length of the density-wave oscillations ( obtained from an 
exponential fit. Arrows mark the corresponding values in the main 
figure. Other parameters chosen as in Fig. [3] 

struction of the incompressible state. This is also signaled 
by an increasing occupation of the A cavities (squares). In 
this regime, the projected model becomes invalid and the full 
numerics very costly as the local Hilbert space cutoff needs 
to be increased substantially. The interesting details of this 
crossover are subject of future work. 

We now investigate the spatial order of the steady 

state by studying the second-order coherence function 

( 2 ) 

(density-density correlator) of the B sites, i.e., g),- n = 

{b\bjbibj)/{b\bi){tfjbj). Fig . jdjshows the spatial correlations 
of the central B site (i = 0) with its neighbors as calculated 
with iTEBD. At weak and intermediate pump strength //ft we 
find local (j = 0) as well as nearest neighbor (j = ±1) anti¬ 
bunching, which represent a signature of photon blockade and 
incompressibility, i.e., the resistance of the system to accept 
simultaneously two pump photons entering the chain either 
on the same or on neighboring plaquettes (which share a qubit, 
see Fig.|T|». Thus, if a photon is present at a B site of the chain, 
every other B site is less occupied due to effective photon- 
photon interactions resulting in polaritonic density-wave like 
order. At larger distances, density-wave order manifests it- 
self in correlations alternating between bunching {g^ 2 j)s > 1) 

and anti-bunching ( 3 o( 2 j+i) B < 1) with a P er i°d doubling of 
two unit cells, leading to an incipient crystalline state of light. 
This can be interpreted as the non-equilibrium counterpart of a 
charge density wave appearing in the ground-state of an elec¬ 
tronic or atomic system with a flat lowest-energy band, e.g., 
in a sawtooth or Kagome lattice ll32ll . 

In a regular one-dimensional Jaynes-Cummings array inter¬ 
actions vanish when g <C J (EL Interestingly, for the flat band 
the converse is true, as the ratio g/J determines the polaritonic 
nature of the plaquette states, which are qubit-like and thus 
strongly interacting when g -C J (see Eq. 0)- This remar- 
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FIG. 5. (color online), (a) On-site and nearest-neighbor correla¬ 
tor SooBiffoiB as a function of g/J for f /n = 0.05 (blue arrow 
in Fig. [3}. The inset shows the complete spatial dependence of the 
coherence function for the g/J values marked by arrows in the main 
figure, (b) Correlation length as a function of g/J obtained from an 
exponential fit of the correlator. Other parameters chosen as in Fig. [3] 


kable effect determines the spatial extent of the density-wave 
correlations. As shown in Fig. [5] the correlation length shrinks 
when the flat band becomes photon-like (g Jp J), while it 
grows steeply in the opposite limit (g <C J). On the tech¬ 
nical level, the projection on the flat band (see SM) modifies 
the drive strength according to f/n — > (//k)\/1 + 2 J 2 /g 2 , 
thus effectively increasing the drive strength when g/J de¬ 
creases. This entails a larger contribution of the high-energy 
density-wave state I’I'dw) (with infinite correlation length) to 
the NESS. At the same time the gap to the other bands closes 
as g/J —> 0 leading to a destruction of the photon blockade. 
Consequently, we find the strongest anti-bunching for the fully 
mixed polaritonic case when g ~ \/'2J. 

Finally, we show that the signatures of geometric frustra¬ 
tion, incompressibility and crystalline order vanish when the 
flat band becomes dispersive, i.e., when <5 QB / 0 (the drive 
stays resonant with the top of the band). For <J QB = J (com¬ 
pare with the dashed lines in Fig.[2ja)), we observe in Fig.[3}a) 
and[4]b) that the plateau as well as density-wave like correla¬ 
tions completely disappear. The latter are replaced by a broad 
and rather featureless anti-bunching dip in Fig. |4jb). Indeed, 
the correlation length of the density-wave oscillations drops 
quickly from its maximum flat band value to roughly one unit 
cell (see inset of Fig. [4]). 

In summary, we have shown that geometric frustration in a 
photonic lattice pushes the system towards an incompressible 
state of light characterized by short-ranged crystalline order 
with period doubling. We have proposed the simplest model 
of a frustrated quasi-ID lattice based on a circuit QED ar¬ 
chitecture realizable with state-of-the-art technology and eas¬ 
ily extensible to two dimensions, e.g., to study topological 
effects. A variant suitable for a realization of our proposal 
based on semiconductor micro-pillar arrays 02 ] [381 is de¬ 
scribed in the supplemental material and has recently been 
realized experimentally 1621 . The onset of long-range corre¬ 
lations motivates another interesting question for future work, 
i.e., whether super-solid behavior of light (coexistence of su¬ 
perfluidity and density-wave order) could be observed in a 
flat band without the need of explicitly engineering nearest- 


neighbor interaction terms in the Hamiltonian fill 1631. Our 
proposal thus paves the way for quantum simulations 1641 of 
frustrated systems far from equilibrium and the realization of 
strongly correlated, exotic states of light with non-trivial spa¬ 
tial correlations. 
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